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A LINEARIZED APPROACH TO CHEMICALLY GENERATED 


WAVES IN A DILUTE, ISOTHERMAL ATMOSPHERE 

PART I 

ANALYTICAL TREATMENT 


INTRODUCTION 

Waves in a stratified fluid under the influence of gravity appear to have 
been initially discussed by Burnside (1889) and Love (1891). Both authors treated 
an incompressible fluid. Gortler (1943) used schlierin photography to show 
experimentally that disturbances in an incompressible stratified medium under 
the influence of gravity propagate along characteristic rays. Lamb (1908) 
treated a compressible, adiabatic, ideal gas whose density is stratified by gravity. 
Since the original contributions by Burnside, Love and Lamb there have been many 
theoretical papers on various types of gravity waves. Reasonably up-to-date 
treatments of the subject may be found in Eckart (1960) and Yih (1965). However, 
gravity wave theory is currently undergoing an active phase of development 
and contributions are being added to the literature at a rapid rate. 

Recent experimental data from the earth's atmosphere and oceans have 
greatly stimulated interest in gravity- acoustic waves. Gossard (1962) has ob- 
served gravity waves in the troposphere. Hines (1960) has shown that internal 
atmospheric gravity waves may account for many of the phenomena observed 
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in the lower ionosphere, and gravity waves have been observed in the thermo- 
sphere by Newton et. al., (1969) and Harris et al., (1969). 

It seems to be generally agreed that gravity waves are generated in the 
troposphere, then propagate upwards. Thus Gossard (1962) has observed 
gravity waves near the earth's surface and Eberstein (1970) has illustrated 
the development of gravity waves between 30 km and IzO km. 

The simpler forms of gravity wave theory treat propagation through an 
inviscid, non-conducting gas, while more sophisticated theories include effects 
of heat transfer, viscosity, and high altitude phenomena such as ion drag. How- 
ever, the effect of chemical reactions does not seem to have been adequately 
considered to date. 

Between the troposphere and the thermosphere there are several regions 
where important chemical reactions take place. 

First, there is the ozonosphere between approximately 15 km and 45 km, 
with peak ozone concentration at about 35 km. Ozone mole fraction is in the 
order of parts per million (Mitra, 1952). 

Between 60 km and 96 km there is Nitric Oxide (Pearce, 1969). The mixing 
ratio for Nitric Oxide is also in ppm. 


Between 90 km and 120 km molecular oxygen dissociates into atomic oxygen. 
The dissociation ratio, a being 3 x 10~ 8 at 90 km and 0.998 at 120 km (Mitra, 
1952). The oxygen is no longer a truly dilute reactant since its mixing ratio is 
20 %. 


The major chemical systems are tied to the absorption of ultra lolet radia- 
tion from the sun s and are thus having a dawn and dusk dependence. There is 
also a dependence on solar activity. 

Large storms and such phenomena as aurorae may also be simulated by a 
chemical type heat release. 

In what follows, an initial study is made of chemically generated waves in a 
dilute, isothermal atmosphere. Simplifying assumptions are made to make the 
mathematics more tractable 
2. Governing Equations 

The continuity and momentum equations (cf. Shere and Bowhill, 1969) are 
not affected by the chemical reaction provided that the state variables are given 
the extended definition: 

P . P (.o, T, a) 

where P, p , T and a denote pressure, density, temperature and degree of 
advancement of the reaction, respectively. These equations are in tensor 


notation: 
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We denote time, space variable (x 2 vertical), velocity component, gravitation 
vector component, viscous stress tensor, dynamic viscosity and the Kronecker 
delta by t, x i , u., g., o\ . , p and 8.. respectively. We also take g A =-gS, 2 
and we define O/Dt to be the total or Eulerian derivative. 

In order to develop the energy equation, it is necessary to consider the 
thermodynamics of the system. The derivation of the equation of state parallels 
the work of Liepmann and Roshko (1957, p. 29) and the derivation of the reaction 
rate equation parallels the work of Eberstein (1966). The following discussion 
applies to the reaction 


J + AB z A+B+J 

where J is an inert element. We assume that each component of the mixture is 
a perfect gas. The mixture, however, is not a perfect gas. The partial pressure 
of the ith constituent is 
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where Rj ' 'R /\V ( with 'R the universal gas constant, W. the molecular weight, 
and rm the mass fraction of ith-component. Letting a be the fraction of AB 
dissociated, n. the number of moles of ith component, and n 0 the total number 
of moles when a = 0, 


AB 

" n o 

X 0 (1 ' 

n A 

O 

C 

i: 

X 0 a 

n B 

n o 

X 0 a 

"J 

= (1 

" X 0 )n 


where X 0 is the mole fraction of AB when a - 0. Summing yields 

= + aX o>* 

P = R(a),oT 
R(a) - JR/M 


(2.4) 


(2.5) 


where !R is the universal gas constant, and M is the mean molecular weight, 
given as: 


M 


( 


W ; + X n (W An - W : 


AB 



1 + a X, 


) 


It follows that 


R(a) = R(0) (l + a X 0 ) 

The reaction equation is as follows: 


Dn AB /Dt 




( 2 . 6 ) 


(2.7) 
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s+% ml rw* 


where 



and k is a rate constant. In determining (2.7) it has been Impllclty assumed that 
all molecules present may act as an Inert element with equal collision effective- 
ness. Substituting (2.4) Into (2.7) yields 


Da 

Dt 


[pk,(T)/M] 



[K(T)p/M] 



( 2 . 8 ) 


where 


k n ( T ) 
K(T) . Ic^tn 


and M is the mean molecular weight 


For chemical equilibrium, the rate equation reduces to 



where use was made of the equality 


K 
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In atmospheric problems the reactant fraction is generally small, and It 
will be assumed that X Q < < 1, It follows that aX Q < < 1, The rate equation 
may thus be simplified to give: 


J>i 

Dt 


['M'O/M] {l - a} 


(2.9) 


The energy equation is now developed in the usual manner (cf. Liepmann 
and Roshko, 1957, p. 185+) . Summing the rate of change of the internal, kinetic 
and potential energy yields 



v 


oc f ~2 / u . u . + / gx 2 j u. n j dA 
= - I Pn . u . dA + I u ; o- . . dA + 

A A 

where n. is the ith component of the unit normal of the closed surface 4>, V is 
the volume of 4> and A the area; e is the internal energy per unit mass and k is 
the coefficient of thermal conductivity. 

By applying Gauss' theorem and substituting the continuity and momentum 
equations into (2.10), we obtain 
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We now use the first law of thermodynamics, h = e + P/7 , and the continuity 
equation to obtain 


Dh „ DP /yu i rl ( <9T\ 

PJJt ’ T5F 4 \ K 17T/ 


( 2 . 12 ) 


In a reacting gas mixture the enthalpy depends on the temperature and the 
degree of reaction, i.e. 

h = h ( T , a ) 

The rate of change of enthalpy per unit mass is then given by: 


Dh 

DF " 



Da 

x qB dF 


(2.13) 


where B is the enthalpy of reaction and C p is heat capacity. For an endothermic 
(exothermic) reaction, B is positive (negative). Both the enthalpy of reaction 
and the heat capacity have a weak temperature dependence, and an even weaker 
pressure dependence. For simplicity both C p and B will be assumed constant 
throughout this paper. 

Eliminating h between (2.12) and (2.13) one obtains the enthalpy equation 
in the form: 
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(2.14) 


p C. 


DT 

Dt 


X 0 B 


Da 

Dt 



To summarize, we have found the system 



( 2 . 1 ) 

( 2 . 2 ) 

(2.5) 

(2.9) 


DT Da 

Dt + X o p B Dt 



( 2 . 11 ) 


The dependent variables are P, p, T, u. , and a; the independent variables 
are t, x .. The quantities k p , K, c p and B are known parameters of the system, 
X 0 is a small (known) parameter. 


For the remainder of this paper we assume that the atmosphere is inviscid 
and nonconducting. Hence, the governing equations are 


Dt 


+ 


<9u 



0 


(2.15) 
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Dt ' Ki ' d: 


(2.16) 


P /> RT 


(2.17) 


St " " (if) [ l ■ °] 


(2.18) 


„ DT t Da DP 
Dt + X o B ^ ; Dt " Dt 


(2.19) 


3. Nondimensionalization of the Governing Equations 

Since mathematical operations are performed on pure numbers it is proper 
to non-dimensionalize the governing equations. 


Some of the non-dimensionalization is almost trivial, as shown immediately 
below: 


R' 

= R/R(0) 

V 

:: C p /R(0) 

p' 

= P/P, 

M' 

= M/M* 

k 9 
k f 

= kp/O^/M*^ 

K' 

- KpyM*co B 
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where R(0) is the gas constant fom o, and the starred reference level is at 
the bottom of the reaction zone. >•< B is a base frequency. 

t ' = ta; B 

X i = X./H 

where H is scale height 

T' = R(0) T/w b 2 H 2 

B 1 = B/6- b 2 H 2 

P* “ P /p* a B 2 H 2 
It is seen that in the primed system we have 

P' “ p'R'T' 

In general, the non-dimensional system of governing equations looks very 
similar to the dimensional system. 

The quantity H 2 is proportional to the energy contained in a wave 
oscillating with the base frequency, and having amplitude equal to a pressure 
scale height. The quantity Ha^ 2 is a measure of the acceleration experienced 
in such a wave. Since such acceleration must be less than that due to gravity 
we would expect g* > 1. 

One may estimate the magnitudes of the non-dimensional parameters. 
Taking a scale height of 6 km and a base period of 5 minutes, one obtains 


o B H ) 2 1.6 x 10 8 cm 2 /sec 2 ~ 3,74 calories/gram 


and 

o g* ~ 3.5 

« B H ~261 cm/sec 2 , giving ^ 

Further, taking B = 34.4 kcal/gmole, M = 48g/gmole, one obtains B’ ~ 200. 
4. Asymptotic Development 

The non-dimensionalized system of equations will be expanded about the 
parameter X 0 , i.e., we expand each dependent variable in a power series of 
the form: 


f(t, X, z) - f (N) (t, x, z) X 0 N (^.1) 

N = 0 

where higher order terms may be dropped as X 0 - 0. This expansion limits 
the validity of our results to dilute fluids, such as the atmosphere. Other ex- 
pansions [cf. Lax, 1957] can be obtained by inserting a small parameter into 
the initial conditions. These expansions may, for example, be used to study 
the generation of gravity waves by tidal waves. 

Assume that the atmosphere is initially stratified, quiescent, and bounded 
below by the earth, i.e., we define an initial-boundary value problem with initial 


conditions: 


T(0, x, z) - T 0 = 
u . ( 0 , x . z ) n 0 ( i “ 1 , 2 ) 

(4.2) 

p( 0, x , z ) n e” * 

<i(0, x, z ) " 0 

and boundary condition: 

U 2 ( t , X , 0 ) = 0 


Using the equation of state to eliminate P from (2.18) - (2.19) and substituting 
(4.1) into the resultant system yields the zero-order solution 


T< 0) (t, x, z) 
Ui (0 > (t, x, z) 
P (0) (t, X, z) 


T = T 
1 o % 


= 0 


e 


(4.3) 


o.(°) (t, x, z) = 1 - exp [-k F/ o( 0 ) t] 


Defining T < 1 ) = T (0) T. 


(i) 
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(i) = o<0) 


p ^ and oS 1 > 


= a<°> 


' a 


( 1 )’ 


we get the first order system 


Bu^ 1 ^ 
Bt 


*P( l 


(i) 


dT 


+ t o-7T^ + 


(i ) 
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(4.4) 


J 2 ^(l) 1 ) ^ Ba< 0) 

at— + t (0) t (1) + T „— 5-J— ^ =- T » — 

+ T 0 a<°> 


^<D au i (1) 

-T^ + ~dT- 
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(9 z 
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(4.5) 


(4.6) 
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_li> + (1) - 

d t + u 2 ( } " 


n ^ (0) 
“® at 


a (D (n) . d«<°> 

at a< > + a ( i ) “at 


- r< n (p<°>, a<°> ) 


where 


ft - ■)• 


Equations (4.4) - (4.7) form a linear system of four partial differential 
equations in four unknown functions. The system is inhomogeneous. The solu- 
tion of (4.8) for a (1 j is needed only for the calculation of second order terms. 
Since we determine only the first order terms in this paper, we will neither 
specify rC 1 ) nor solve for a . 

Dropping the (1) subscript or superscript and eliminating we obtain the 
equation 


Ill 

'a 2 

- c 2 - c 2 ^4- 

+ c.il 


[at 2 

_at 2 

d \ 2 az 2 

azl 

r 2 ax 2 J 


f(a<°>) = - (y- 1) 13 


Taj 
t 2 [at 


! a 2 a 2 

1 ~ T ° aT 2 " T ° dP a(0> 


a 3 a(°) 

<r " 1 )B T^ 7 T + 
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a 2 a<°> 




and 


> T o 


After (4.9) is solved, the other dependent variables can be determined by 
the equations 


Lf. * *] 


_ , , 1 c) 2 T ,7T „(9 2 a<°> . ,, ,,, 

+ t„ x; I p = j—i-J-f - T oT7 + T 0 T + i„a<°> (4.11) 


<9u 


u i _ d l p i d 2r r i aT 

"i *■ ■* — j i * “* r ” .. «■> 


dx dtdz y ~ l dt dz y 




and 


<9/0 1 <9T aa<° ) 

u 2 ’ TT " ITt " 8 Tt 


(4.13) 


Now define 


^(t, x, z) = e~ z/2 T(t, x, z) 


(4.14) 


The above definition allows the temperature perturbation to grow expo- 


nentially with altitude without a corresponding growth in 6 . 


Substituting the definition of 6 into (4.9) yields: 


a 2 

at 2 



c 2 v 2 



c 2 (y i) a 2 e 


>-*/2 


7“ 


ax' 


f(a) 


(4.15) 


f(a (0) ) was defined in Equation (4.10). 


#w 
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The dispersion equation of (4.15) is: 


x t 4 " C 2 ( A « 2 + ^. 2 1 t) \ 2 + [(■>" !) CVy] \„ 2 = 0 (4.16) 

where is the wave number of i. 

For a discussion of dispersion relations and how these are obtained from 
partial differential equations reference is made to Courant and Hilbert, es- 
pecially p. 588. 

5. Solution of the First-Order System 

From (4.4) -(4.6) and the initial conditions (4.2) we deduce that 


0(0, x, z) = 0 

(5.1a) 

0 t ( 0, x, z) = - (y - 1) Bk f e' 3 * /2 

(5.1b) 


Conditions for 0 tt and 6 ttt at t = 0 can be determined by differentiating (4.4)— 
(4.8) appropriately and analyzing the resultant system or by integrating (4.15) 
and applying the physical requirement that <9(t, x, z) must remain bounded as 
t - oo. The negativeness of & t (0, x, z) means that the atmosphere initially 
cools, which is in correspondence with the endothermic reaction chosen. 

Since the inhomogeneous parts of (4.5)- (5.1) do not depend upon the hori- 
zontal space variable, neither does the solution. Thus (4.15) reduces to 
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(5.2) 

(5.3) 

It is worth noting that the initial conditions and the alpha dependence chosen 
are such as to imply the relevance of only one space dimension. Ot tc* initial 
or boundary conditions would require that the second space dimension be kept. 
Integrating (5.2) with respect to time twice, 

M [ fl ] = w(t, z) = k F (y - 1) (k f [- R + Bt 2 ] e -z - 2Bt) exp - 

subject to initial conditions (5.1). 

The inhomogeneous part of (5.4) is the effect of the chemical reaction. For 
a realistic approximation to the atmosphere, the inhomogeneous part of (5.4) must 
be multiplied by a suitable weighting function, since the chemical contribution 
only extends over a finite altitude regime. We also note that at each fixed 

4 

altitude the driving force must tend asymptotically to zero in the Poincare 
sense as t - oo . 

The assumptions used in this paper are also suitable for a multilayer ap- 
proximation of the atmosphere. In this case we would consider either an initial- 
boundary or pure boundary value problem for (4.15). The x-dependence may not 
generally be eliminated for this case. 


3z 

T 


+ k F te’ 


(5.4) 
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The solution of (5.1)-(5.4) is obtained by superposition. We choose a suitable 
function tt (t, z) which satisfies the initial conditions and is damped in time. 

For example, set 

'^(t, z) r ~ - (7 - 1) 8 k f te~ v/2tc e ” 3t/2 (5.5) 

Observe that (5.1) is satisfied und 

M [ 1 ,^) - 2/2 C (7 - 1) 3k f e _/2 c ~ 3 * /2 

Since M [ ] is a linear operator, (5.1)-(5.4) may be reformulated as 

0(t, z) - 0 f (t, z) + fJ 2 ( t, z) (5.6) 

where () 2 (t, z) satisfies 

M[0 2 (t, z )] w(t, z) - 2/2 C (7 - 1 ) Bk f e' /2 tc e" 3l/2 , (5.7a) 

and 

0 2 (0, z) n 6 ^( 0 , z) = 0 . (5.7b) 

The solution of (5.7) for 0 2 ( t, z) is known Ccf. Lamb, 1909] and is given by 

0 2 (t, z) = f W(t, z; r) dr (5.8) 

Jo 

where 


/* t C"T C I 

Jo (I 1 

V “ 

rf 

N 

HI 

/(tc - re) 2 - 77 2 J 

^0 1 


/ 


•Q(r, z, 77 ) dTy (5.9) 
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and 


with 


QO, z, v) I . ( \ z ( •'/) i ?0r , z - '/)] ''2c 2 
K*. z) w(7, z) - 2 t / 2c( / - 1) Rk f (’ v ' 2 ' tf c" 3l/2 


(5.10) 


Although (5.5) and (5.8)-(5.10) provide an exact solution of (5.1)-(5.4), this solu- 
tion is complicated and yields quantitative information only through cl merical 
calculations. The results of these calculations are given in the following sections. 


A solution of (5.1)-(5.4) can also be obtained using series techniques. The 
result is a somewhat rapidly convergent infinite series whose terms involve 
products of powers (- k f te"*) n and Young's functions [Young, 1912], a special 
case of Lommel functions. Since the series was not used to obtain physical 
results, it is not presented here. 

From (4.12) -(4.13) we deduce the first order terms: 

T(t, z) = e z/2 0( t, z) 


r z $ ~ t 

U 2 (t,z) = - I jj: + dz 

J - Z L 


where z is at the earth surface and 

s 


p{ t, z) = 


f [- - *] 


dt 


The numerical solution and appropriate physical applications are discussed in 


Part II. 
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